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Abstract 

We address the question about the velocity of signals carried by 
Bessel beams wave packets propagating in vacuum and having well 
defined wavefronts in time. We find that this problem is analogous 
to that of propagation of usual plane wave packets within dispersive 
media and conclude that the signal velocity can not be superluminal. 
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1 Introduction 



It is generally accepted that the causality principle in special relativity im- 
poses the speed c of light in vacuum as the upper limit for the velocity of 
propagation of signals or interactions. By the way, Maxwell equations in ma- 
terial media allows, in special circumstances such as anomalous dispersion, 
plane wave solutions propagating with superluminal (faster than c) phase 
velocities M . Although these velocities cannot be readily identified with the 
signal velocity, the question about the possibility of building a superluminal 
signal by superimposing these solutions naturally arises. 

The above problem was studied by Sommerfeld and Brillouin at the be- 
ginning of the last century ||. In order to define precisely a velocity of prop- 
agation Sommerfeld considered only signals carried by wave packets which 
have well defined wavefronts (and possibly ends) in time, which we call SB 
signals. The signal itself, or its "main part", independently of the way it 
is defined, is obviously confined to the region behind the wavefront. These 
packets were assumed to be normally incident on a semi-infinite dispersive 
medium, in which boundary the arrivals of both the wavefront and the signal 
were assumed to coincide. The Sommerfeld's main result was that, irrespec- 
tive the material medium, the wavefront of such a packet propagates always 
with the speed c. Consequently, the arrival of the signal at a given point in- 
side the medium can occur only after that (or simultaneously, if the medium 
is nondispersive) this point is reached by the wavefront. Therefore, Sommer- 
feld result implies that the speed c of light in vacuum is an upper bound to 
the velocity of propagation of SB signals. Brillouin studied in great detail 
the evolution of these wave packets within dispersive media, giving a precise 
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definition for the signal velocity, which agrees with the Sommerfeld result. In 
what concerns the group velocity, which can be greater than c in the presence 
of anomalous dispersion, or even negative, Brillouin stated that "...the group 
velocity has a meaning only so long as it agrees with the signal velocity." 

Recently, there is a growing interest on the subject of superluminal wave 
motion in the literature. For example, we can quote experiences measur- 
ing superluminal velocities in the passage of wave packets through barriers 
(see, for example, ||-|| and references therein) and a considerable amount 
of works, both in theoretical and experimental views, concerning the super- 
luminality of Bessel beams and X waves (P|-|[15|| and references therein). 
These works raised questions about the interpretation of such superluminal 
velocities, specially in what concerns the meaning of signal velocity and its 
connection with the group velocity and the causality principle. So, this is 
still a very debated subject. 

In this letter we are concerned to Bessel beams wave packets, from which 
X waves are a special case. These beams are (inhomogeneous) plane wave 
solutions of the homogeneous scalar wave equation in vacuum which propa- 
gate with superluminal phase velocities. X waves are localized waves built 
up as special superpositions of Bessel beams and propagate rigidly (without 
dispersion) in vacuum with superluminal velocities. From these nondisper- 
sive properties some authors suggested that this superluminality could also 
be associated with the signal velocity |L0], [11], We analyze the question 
of the maximum velocity of signals carried by Bessel beams wave packets. 
To this aim we follow the approach of Sommerfeld and Brillouin cited above 
and consider only SB signals, i.e., those carried by Bessel beams wave packets 
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having well denned wavefronts (and ends) in time. We first study the chopped 
Bessel beam, which have a finite duration at its source. From the analysis 
of this packet we identify a mathematical analogy between its propagation 
properties and those of usual (inhomogeneous) plane wave packets propagat- 
ing within a dispersive medium, namely a tenuous electronic plasma. This 
is our main result, because such analogy makes possible a straightforward 
application of Sommerfeld result to conclude that the wavefronts of these 
packets propagate with the speed c, while the wave packets distort while prop- 
agating. This is a curious result because we are dealing with propagation of 
waves in vacuum. As a direct consequence, the velocity of the signals carried 
by these wave packets can never be superluminal, independently of the way 
it is defined. 

We also consider briefly the experiment of Mugnai, Ranfagni and Ruggeri 



Tl| , which posed a question about the superluminality of signals carried by 



Bessel beams. An explanation for the measured velocities in this experiment 
was given in references [|T3| , p~5f in terms of interference phenomena, showing 
that the superluminal velocities of the peaks moving along z axis were not 
causally connected and, therefore, did not represent signal velocities. To be 
able to say something about the signal velocity in this experiment we argue 
that the waves produced in it can be viewed as a kind of finite aperture 
approximation to chopped X waves. The last ones are ideal waves (they 
need an infinite aperture to be produced) built up as linear superpositions 
of chopped Bessel beams. We show that chopped X waves have wavefronts 
moving with velocity c and thus the SB signals carried by them can not have 
superluminal velocities. Also, we suggest that the observed superluminal 
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peaks can be qualitatively explained from the fact that the chopped X wave 
distorts along the propagation, showing a kind of reshaping phenomenon, as 



observed in 13 



The generalization of the analogy with dispersive media to a large class 
of superpositions of Bessel beams is straightforward and will be done in 
Section 4. In the last section we make our concluding remarks. In particular, 
we comment on some discrepancies between our results and others in the 
literature concerning the wavefront velocities of chopped Bessel beams and 
chopped X waves. 

2 Bessel beams wave packets 

Bessel beams are cylindrically symmetric solutions of the scalar homogeneous 
wave equation in vacuum. They are given by @, §, pi], 



¥ k ^\p,z,t) = J (pfc p )exp{«(fc,*-wt)}, (1) 

where Jq(x) is the Bessel function of the first kind and order zero and the 
parameters k z , k p and u satisfy the following relation 



— o 
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uj — c k„ 



(2) 



Thus, from a mathematical point of view, any two parameters among k z , k p , 
and u) can be chosen independently. For the purposes of this letter we choose 
k p and u as the independent ones and assume they are real. Making so, k z 
will be given by (§) and it can, in principle, be imaginary. 

At this point we emphasize that the above relation is mathematically 
identical to the dispersion relation of a tenuous electronic plasma if we iden- 
tify c 2 k 2 p with the square of the plasma frequency u 2 [|TJ. 



If uj 2 > c 2 k 2 p the Bessel beam represents an unidimensional wave 
motion propagating along the z direction, whose propagation properties are 
governed by the real wave number k z . The surfaces of constant phase are 
planes perpendicular to the z axis which propagate with the phase velocity 
v p , given by 

Vp = F ' ^ 

such that \v p \ > c, i.e., the phase velocity in vacuum is superluminal. These 
waves are also inhomogeneous plane waves, due to the presence of Bessel 
function Jo(pk p ) in ([[]). In the case k p = the Bessel beams degenerate 
to usual homogeneous plane waves, with \v p \ = c. On the other hand, if 
uj 2 < c 2 k 2 the wave number k z is imaginary. In this case there is no wave 
propagation in the z direction, but instead of a behavior analogous to ab- 
sorption or attenuation. 

The most general (complex) solution formed from superposition of Bessel 
beams ([!]) is given by 

z, t) = V + (p, z, t) + z, t) , (4) 

where 

roo roc 

ty ± (p,z,t) = / dk p duA ± (kp,u)J (pk p )exp{i[±k z (k p ,Lj)z—ujt]}. (5) 

JO J — oo 

In this expression A^^kp^u) are spectral distributions and the dispersion 
relation is given by [|J 



if U > ck n 



k z (kp,u) = { <• v , " 2 ' "^i-^ . (6) 

- ' ' M.f^KZi, if u <ckn 



With this choice of signals the terms ^ + (p, z, t) and *&~(p, z, t) in the super- 
position (ffl) correspond, respectively, to right and left moving wave packets. 



Specifying the boundary conditions at z = 0, 

*(p,o,t) = *°(p,t); (7) 

— (p,0,t) = *'°(p,t), (8) 
we can determine the spectral distributions A ± (A;p,a;) in (|])[] 

/"OO /"CO f 1/1 

A ± (k p ,u) = -^ dte^ dppJ (pk p )U°(p,t)± s 9°(p,t)\. 

2tt J-oo Jo { ik z {kp,u>) J 

(9) 

To end this section we consider the situation in which only propagating 
modes of Bessel beams (P enter into superposition i.e., we restrict ck p < 
\u\ in (|5T). In this case we can write 



\lu\ u 
k p = — sin 9 and k z = — sin 9 , with 9 G 

c c 



0,|). (10) 



Using these relations we can change the integration variables in (|]) from 
{fc p ,u;} to {9, ZJ} according to 



uj = lo ; 



p 



9 = arcsin 

In this way the two terms in superposition (|4]) can now be written asQ 

^(AM) = [* d0 [°° dujB ± (9,u)J (p-sm9) e i ^^ zcosd - ct \ (11) 

JO J-oo \ C J 

1 To derive this result we make use of orthogonality of Bessel functions 

f°° 1 
/ dppJ (pk p )J (pk p ) = TrS(k p - k' p ) . 
Jo K p 

2 We are taking into account the fact that J is an even function. 
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where B ± {9^uo) are the spectral distributions with respect to the new param- 
eters. The parameter 9 is called the axicon angle. Clearly the superposition 
is more general and includes this last one. 
When the axicon angle 9 is fixed to some value 9 , the last expression 
yields 

*l(j>, z, t) = duJS^Jo [p- c sinfl ) e m±z^o-ct] (12) 

This expression is frequently found in the literature and defines the so-called 
X waves.0 Superposition ( |T2"D can be viewed as a cylindrically symmetric 
superposition of plane wave packets tilted over the z axis by an angle 6q || |[| . 
The particular form of these packets is given by the spectral distributions 
S ± {uj). As an example, if we choose 

=1 , (13) 
to . 

then the plane wave packets are given by a rectangular pulse modulation of 
a carrier of frequency u c . In this expression T is the time duration of the 
pulses and to is the instant of time in which their wavefronts reach the origin. 

At this point we observe that X waves ( |T2] ) are infinitely extended along 
the direction of propagation z. So, they do not have a well defined wavefront 
into this direction and thus do not define SB signals. In the next section 
we shall consider the superposition ([§-[5]) which "chops" (perpendicularly 
to z direction) a Bessel beam characterized by a given frequency cuo and 
axicon angle 9 . From these chopped Bessel beams we will construct the 
corresponding chopped X waves. 

3 If S + (uJ) = S(uJ — l)q) then defines a Bessel beam (|l|) characterized by the axicon 
angle 8q and frequency luq propagating to the right. The same occurs if S~ = 5, with 
motion to the left. 



S ± (u) 



-ic 



An (to 



UJr 



j(u)—iOc)To 



- 1 
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3 Chopped Bessel beams 



We now consider a SB signal consisting of a Bessel beam modulated in time 
by a rectangular pulse at the plane z = 0, i.e., a chopped Bessel beam. This 
signal is introduced in ]T3| through the boundary conditions ([7|-§]), with 



¥°(p,t) = T(t)J (p^8mflo)e-*** (14) 



*'°(p,t) = T(t)(^cos0 o J Jo^sin^oje-^, (15) 

where o;o and #o are /ize<i constants and the time modulation is given by 
T(i) = O (£) — ©(£ — r), with 0(£) being the Heaviside step function. To 
determine the spectral distributions in the general superposition (§-^D which 
correspond to these boundary conditions we substitute these last expressions 
into (||). Then we have 

A ± (k p ,u) = 5^ p -^sin^o)c ± H, (16) 

where 

1 e^-lf : fcosflp | 

We see that the spectral distributions A ± (k p , uj) split into a factor depending 
only on k p and another factor depending only on u>. Then, we can write the 
right moving term in (^) as 

+ (p, z, t) = J™ dk p 5 (k p - ^ sin O ) </>+(A; p , p, z, t) , (18) 

where 

/oo 
duj C + (uj)J (pkp) exp{i[k z (k p , uj)z — ut]} . (19) 
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We can readily identify this last integral with that describing a wave packet 
formed from superposition of (inhomogeneous) usual plane waves and moving 
within a tenuous electronic plasma, with plasma frequency given by uj p = ck p . 
As the spectral distribution C + (uo) is analytical in the upper half of the 
complex cj-plane, this packet itself has a well defined wavefront in time at 
z = [0, @. Therefore, we can readily apply the results of Sommerfeld and 
Brillouin theory to conclude that this wavefront moves with the velocity c 
and the wave packet distorts while propagating. 

To completely describe the wave motion to the right, we must now con- 
sider the k p integral in ([TJ). So, we have 

*^p,z,t) = <j> + (^-,p,z,tj , (20) 

where in this expression to® = u sin9 . Thus, all the conclusions after equa- 
tion fll^) remains valid substituting the plasma frequency ck p by cu®. These 



conclusions also hold for the left moving part of the wave packet (|J) . 

Superposing now the chopped Bessel beams above via spectral distribu- 
tions S ± (uq), we have 

< op (p, z, t) = *t )chop (p, z, t) + v^ )chop (p, z, t) , (21) 

where 

¥ e f chop (p,z,t) = [°° duoS^Uo)^ (^p,z,t) . (22) 

It is easy to verify that in the absence of the time modulation T(t) this 
superposition would satisfy the same boundary conditions as the X wave (p~2|) 
if we identify the spectral distributions in the two expressions. Accordingly, 
we call superposition (|2l| - |2"2"D a chopped X wave. As this wave is a linear 
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superposition of chopped Bessel beams, which have well defined wavefronts 
propagating with the velocity c, also the chopped X waves will have this 
property. As a consequence, they cannot carry superluminal SB signals. 

Now we consider briefly the experiment of Mugnai, Ranfagni, and Rug- 
geri, in which the authors measured superluminal velocities in the propa- 



gation of X waves produced experimentally by a finite aperture device [11 



They posed the question about the possibility of interpretation of these su- 
perluminal velocities as velocities of signals. Causal explanations for the 
results of this experiment were given in references fll3| [15| based on simple 
models showing interference phenomena of waves produced outside the axis 
along which the superluminal velocities were measured. These explanations 
refute the possibility of interpretating the measured superluminal velocities 
as velocities of signals in that experiment. In the reference [13[] the authors 
call the attention to the fact that this experiment shows a kind of generalized 
reshaping phenomenon occurring in free space, characterized by the fact that 
the peak travels along the symmetry axis faster than the wavefront. 

In order to apply our analysis to this experiment we assume that each 
pulse of the experimentally produced waves can be viewed as a finite aperture 
approximation to an ideal chopped X wave. We consider the propagation 
in the region z > 0, where the plane z = contains the borders of the 
mirror which produce the wave during the time interval [0, r]. To set the 
theoretical (infinite aperture) model to the experiment we consider only the 
right moving packet in (12), with the spectral distribution S + given by QI^I). 
Now we interpret the parameter u c as the frequency of the microwave carrier 
in the experiment, T as the time duration of the rectangular modulation of 
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this carrier and t = Rs ™ e ° as the instant in which the peak (which travels 
along the z axis) begins to be generated at the z = plane. R is the radius of 
the aperture and 9q is the axicon angle fixed by the experiment .f\ In order to 
produce the appropriate boundary conditions in the region of finite aperture 
(p < R) at the plane z = 0, the time duration r of the chopped X wave must 
satisfy r > T + 2t . 

Assuming this model, our results about propagation of a chopped X wave 
imply that its wavefront propagate with velocity c and can not carry a su- 



perluminal signal, a result that agrees with references [JTS], [15fl . Also, we 
suggest that the cited reshaping phenomenon could be at least qualitatively 
explained from the fact that the wave packet distorts along its propagation. 

4 General results 

In the last section the essential feature allowing us to identify the analogy 
between the propagation of the chopped Bessel beam in vacuum and the 
propagation of usual plane wave packets within a dispersive medium was 
the factorization of the spectral distribution A (k p ,u)) into a product of a 
distribution depending only on k p and another depending only on to. We can 
generalize our results for the class of all spectral distributions satisfying this 
property. For simplicity, we will concern us to the right moving part of the 
wave packet (f|). All the conclusions will also hold for the left moving part. 
Let A + {k p ,uj) = A p (kp)A w (u). Then, from fl5|) we have 

poo 

+(p, z,t)= dk p A p {k p )<p + {k p , p, z, t) , (23) 

•J o 



4 The axicon angle is given by 9q = arctan jj , where d is the mean diameter of the slit 



and / is the focal length of the mirror. 
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where 

/oo 
du Au(u) exp{i[k z (k p , u)z - ut]} . (24) 
-oo 

Following the same lines of the last section, we identify this expression as 
representing an inhomogeneous usual plane wave packet propagating within a 
tenuous electronic plasma, with plasma frequency ck p . Again, the k p integral 
in (|23|) denotes a superposition of these packets. If each component packet 
have a well defined wavefront in time at some plane perpendicular to the z 
axis, the complete superposition will also have a time wavefront which, by 
our previous analysis, propagates with velocity c. 

To summarize, the velocity c is the upper bound for the velocity of any 
SB signal carried by a wave packet belonging to this class. 

5 Concluding remarks 

In this letter we studied the wavefront velocity of time limited (chopped) 
Bessel beams wave packets and showed that they propagate in vacuum in a 
way analogous to usual (inhomogeneous) plane waves propagating within a 
tenuous electronic plasma. From this analogy we were able to apply the Som- 
merfeld and Brillouin results, originally conceiving the propagation of usual 
plane wave packets within dispersive media, to conclude that the wavefronts 
of these wave packets move always with velocity c, while the waveform dis- 
torts along the propagation. These conclusions were generalized for a large 
class of wave packets having well defined wavefronts and described by fac- 
torizable spectral distributions. 

The above results contradict some conclusions of references fO, fl3, in 
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which the authors conclude that chopped Bessel beams propagate without 
distortion with superluminal velocity. By a careful computation of the spec- 
tral distributions in (H) which correctly give the boundary conditions ( |l|j - 15) 
we found that the spectral distribution for the parameter k p is a delta dis- 
tribution which fixes this parameter at the value ^sin6 , contrary to the 
claims of the authors of |13| , which consider k p as a parameter depending 
on the varying frequency uj (not to be confused with the fixed parameter co> , 
characterizing the boundary conditions)^. 

From superposition of chopped Bessel beams we constructed chopped X 
waves and suggested that these waves could be viewed as infinite aperture 
theoretical models for the waves produced in the experiment of Mugnai et al, 
which posed a question about the superluminality of signals carried by Bessel 
beams. By our analysis the chopped X waves have wavefronts which move 
with velocity c. Therefore, they can not carry superluminal signals. Also, 
from the fact that these packets distort along the propagation, we suggest 
that the reshaping phenomenon cited in Jl3| could be at least qualitatively 
explained. It would be interesting to develop the analogy with dispersive 
media further in order to compare also quantitatively the results arising from 
this model with the experimental data of Mugnai et al and with the models 
based on spherical waves presented in (13], |15| . 

Our analysis introduced the formal analogy with dispersive media as an 
alternative way to approach the problem of velocities of signals carried by 
Bessel beams. We hope it can help us to better understand this debated 
question. 

k p , in our notation, corresponds to the separation constant fl, in their notation. 
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